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FOREWARD 


The  spectral  sea/swell  analysis  and  forecast  model  described 
in  this  technical  note  is  the  first  to  be  adapted  for  real-time  use. 

The  program  has  been  running  twice  daily  in  the  FNWC  operational 
program  stream  since  April  1972.  It  was  developed  through  the 
cooperative  efforts  of  the  Naval  Oceanographic  Office,  New  York 
University,  and  Fleet  Numerical  Weather  Central. 

In  working  with  any  spectral  model,  verification  is  always  a 
major  problem.  Spectral  data  are  quite  scarce;  calibration  and  tuning, 
etc. ,  are  very  difficult.  Fortunately,  in  this  case,  we  had  the 
services  of  the  Naval  Oceanographic  Office  representative  to  FNWC, 
who  arranged  the  joint  verification  project  reported  here. 

Future  plans  for  this  model  include  evaluation  of  Environmental 
Prediction  Research  Facility's  (EPRF)  Fields  by  Information  Blending 
(FIB)  wind  analysis  of  the  Mediterranean  "window"  as  input,  as  soon 
as  available. 


Reviewed  and  approved  1  March  1973. 


W.  S.  HOUSTON,  Jr/ 

Captain,  U.  S.Navy 

Commanding  Officer 

Fleet  Numerical  Weather  Central 
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Abstract 


An  operational  Mediterranean  wave  spectral  model  has  been 
developed  through  the  cooperative  efforts  of  the  U.S.  Naval 
Oceanographic  Office  (NAVOCEANO),  Fleet  Numerical  Weather 
Central  (FNWC)  and  New  York  University  (NYU).  Since 
April  1972,  the  model  has  been  operating  in  a  real-time  environ 
ment  (analyses/prognoses  to  48  hours).  The  model  is  a  modified 
version  of  the  original  NYU  North  Atlantic  wave  model  and  has 
two  main  parts:  (1)  wave  energy  growth  based  on  a  modified 
version  of  the  Miles-Phillips  growth  mechanism  and  dissipa¬ 
tion  at  individual  grid  points  and  (2)  wave  energy  propagation 
from  grid  point  to  grid  point.  The  Mediterranean  wave  model 
uses  a  conic  conformal  grid,  permitting  the  assumption  of 
equal  spacing  between  grid  points.  There  are  455  sea  points 
with  a  mesh  length  of  67  km. 

The  wave  model  driving  force  is  a  modified  version  of  the 
FNWC  Marine  Wind  model  which  has  a  mesh  length  of  approximately 
370  kms.  For  use  in  the  Mediterranean,  the  winds  are  inter¬ 
polated  between  Marine  Wind  field  grid  points.  At  analysis  (or 
pre-analysis) times ,  wind  reports  from  synoptic  ship  files 
are  reanalyzed  in  the  Mediterranean  wind  program  so  that 
local  wind  phenomena,  such  as  "mistrals,"  are  included  in 
the  wind  field. 

The  Mediterranean  wind  and  wave  model  was  evaluated 
with  ship  reports,  one  wave  buoy  and  an  airborne  wave  laser. 
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The  wave  computations  appear  to  be  as  good  as  the  wind  input 
allows. 

The  on-time  analyses  and  24-hour  forecasts  of  wave  heights 
are  transmitted  to  Fleet  Weather  Central,  Rota,  Spain. 
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1.  Introduction 

A  Mediterranean  (Med)  Sea  wave  spectral  forecasting  model 
has  been  developed  through  the  cooperative  efforts  of  the 
U.S.  Naval  Oceanographic  Office  (NAVOCEANO) ,  Fleet  Numerical 
Weather  Central  (FNWC)  and  New  York  University  (NYU) .  Although 
various  wave  spectral  models  have  been  in  existence  for  several 
years,  none  of  the  models  have  been  used  for  wave  forecasting 
on  an  operational  basis  primarily  due  to  computer  memory  and 
time  requirements.  Environmental  facilities  were  constrained 
to  developmental  spectral  modeling  and  hindcasting  projects. 
Thus,  the  Med  wave  model  is  the  first  operational  wave  spectral 
forecasting  model. 

The  Med  wave  model  is  a  direct  descendant  of  the  original 
NYU  North  Atlantic  model  (Baer,  1962),  which  has  been  used 
several  times  for  hindcasting  projects,  including  a  six-week 
project  conducted  by  the  NAVOCEANO  representative  at  FNWC 
during  the  summer  of  1971.  At  the  completion  of  this  NAVOCEANO 
study,  it  was  apparent  that  the  North  Atlantic  model  could  be 
adapted  as  a  wave  forecasting  model.  The  only  serious  limita¬ 
tion  of  this  model  was  its  grid  system,  a  Lambert  Conformal 
Grid  with  a  mesh  length  of  200  kilometers.  Since  the  Lambert 
Conformal  Grid  causes  large  distortions  over  great  distances, 
it  was  decided  to  limit  the  operational  use  of  this  model  to 
small  geographical  areas  such  as  the  Mediterranean  and  South 
China  Sea. 
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2.  Wave  Theory 


There  are  two  parts  to  the  basic  wave  program:  (a)  growth 
at  each  grid  point  and  (b)  propagation  from  grid  point  to  grid 
point.  The  wave  model  treats  these  two  factors  independently, 
a.  Wave  Growth 

The  scheme  for  computing  the  growth  of  the  wave  energy 
is  based  on  the  work  of  Inoue  (1967).  Inoue,  under  the  guidance 
of  Professor  Willard  J.  Pierson,  combined  and  modified  the  Miles 
and  Phillips  growth  mechanisms  for  practical  application  to 
wave  forecasting.  Although  Inoue  has  described  this  mechanism 
in  great  detail,  the  highlights  of  his  work  as  related  to  the 
Mediterranean  Sea  will  be  briefly  described  in  the  following 
paragraphs . 

The  Miles  instability  theory  (1957)  stated  that  the 
mean  rate  of  energy  transferred  from  the  parallel  shear  flow 
to  the  surface  wave  is  proportional  to  the  curvature  of  the  wind 
profile  and  inversely  proportional  to  the  slope  of  the  wind  pro¬ 
file  at  the  height  where  the  mean  wind  velocity  is  the  same  as 
the  phase  speed  of  the  wave  component.  This  theory  was  original¬ 
ly  proposed  to  overcome  the  sheltering  coefficient  difficulties 
encountered  by  Jeffries  in  1926.  Jeffries  introduced  the 
sheltering  coefficient  to  improve  the  Helmholtz  instability 
theory  (as  described  by  Lamb,  1932)  which  required  stronger 
winds  than  actually  observed  to  generate  waves.  The  sheltering 
coefficient  tried  to  describe  empirically  the  pressure  differ¬ 
ence  distribution  that  occurs  when  a  wind  passes  over  a  wave, 
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i.e.,  high  pressure  on  the  windward  side  and  lower  pressure 
on  the  leeward  side.  Jeffries'  results  were  an  order  of 
magnitude  too  large.  One  problem  with  the  sheltering  concept 
is  that  it  will  not  work  in  calm  seas  and,  thus,  must  assume 
the  presence  of  existing  waves. 

The  Phillips  theory  essentially  stated  that  a  reso¬ 
nance  between  the  air- sea  system  could  occur  when  a  component 
of  the  surface  pressure  distribution  moved  at  the  same  speed 
as  a  free  surface  wave  of  the  same  wave  number  (where  the  wave 
number,  k,  is  equal  to  2ir/L,  and  L  is  equal  to  wave  length). 
Unlike  prior  growth  theories,  this  theory  can  be  applied  to  a 
calm  sea  surface.  Miles  developed  a  combined  resonance- 
instability  model  in  1960.  Snyder  and  Cox  (1966),  using  a 
wave  recorder  towed  by  a  boat,  and  Barnett  and  Wilkerson  (1967), 
using  an  airborne  wave  recorder,  showed  that  the  wind  waves 
were  affected  by  both  mechanisms.  Inoue  then  combined  the 
work  of  Miles  and  the  more  recent  efforts  of  Phillips,  (1966). 
One  of  the  important  results  of  this  work  is  that  when  the 
sea  starts  growing  from  a  calm  condition,  the  resonance 
mechanism  predominates  and  later,  after  the  wind  has  been 
blowing  for  a  given  time,  the  instability  mechanism  takes  over. 
The  resonance  mechanism  leads  to  a  linear  growth  as  a  function 
of  time  while  the  instability  mechanism  leads  to  exponential 
growth,  as  shown  in  Figure  1.  Cardone  (1969)  claims  that  the 
resonance  term  contributes  only  a  small  part  of  the  energy  to 
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the  waves  and  its  chief  function  is  to  trigger  growth  by  the 
instability  mechanism. 

Neglecting  non-linear  effects,  the  spectral  components 
can  be  initially  expressed  as 


S(f,t,x)  =  A[f ,u(t , x) ]  +  B(f ,u(t,x) ]  •  S(f,t,x)  (1) 

where  A[f, u(t,x)3  is  the  resonance  growth  mechanism, 

B[f,u(t,x)3  •  S(f,t,x)  is  the  Miles  instability  term, 
S(f,t,x)  is  the  spectral  density,  and 

f  is  the  wave  frequency,  t  is  time,  x  is  wave  direction. 
Inoue  showed  that  A(f,u)  could  be  expressed,  for  the  North 
Atlantic,  as 


_ _ A* (to) 5,25  U2-2* _  dQ 

-tt/2  [1/4  (jj)2  +  (k  sin  0)2][l/9  (Jf) 2 * S  +  (k  cos  e  -  g)Z] 


(2) 


where  w  *  2trf, 

U  =  the  wind  speed  at  a  certain  anemometer  height  (NYU 

uses  19.5  meters  as  this  is  the  average  anemometer 

height  of  the  weather  ships), 

2 

k  *  Co  /g  and  g  is  the  gravitational  acceleration, 

A*  is  directly  proportional  to  k  and  inversely  propor¬ 
tional  to  the  water  density  and  the  square  of  gravity. 
A*  was  at  first  determined  from  observational  data  and  then 
modified  as  the  original  NYU  wave  programs  were  used  for  hind¬ 
casting.  A*  is  a  combination  of  all  the  constants  as  well  as 
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some  empirical  calculations.  Inoue  used  the  value  of  9.84 
*  ID'15  for  A*.  This  value  was  calculated  from  the  Snyder 
and  Cox  experiment.  Inoue  tried  to  detect  the  A  growth  term 
from  spectra  calculated  from  British  weather  ships  and  Argus 
Island  Tower  data  but  was  unable  to  reach  any  definite  con¬ 
clusions.  He  suggested  that  the  value  shown  in  equation  (2) 
be  used  and  that  this  growth  value  be  subtracted  from  obser¬ 
vational  data  to  determine  the  (B»S)  term. 

Before  describing  the  B  term,  several  other  relation¬ 
ships  need  to  be  introduced.  The  first  relationship  is  one 
which  describes  the  wind  profile  above  the  sea  surface  and  is 
expressed  by 


U(Z)  -  C/cos  8  =  U*/kQ  In  (Z/Zm)  (3) 

where  U*  is  the  frictional  velocity  equal  to  AQ/pa  and 

t  is  the  surface  stress  and 
o 

p  is  the  atmospheric  pressure, 

«& 

k  is  a  von  Karman  constant, 
o 

C  is  the  phase  velocity 

Zm  is  the  matched  layer  where  the  wind  speed  and 
phase  speed  are  the  same  [U(Zm)  cos  8  *  C) . 

By  combining  several  terms  in  equation  (3),  a  rough¬ 
ness  parameter  ZQ  can  be  expressed  by 

Zo  ■  Am  eItP  (U»  cos  61  (4) 
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Several  empirical  relationships  have  been  derived 
which  describe  the  relationship  between  ZQ  and  U*.  Inoue 
used  the  expression 

Z0  *  U*/g  (5) 

which  is  based  on  the  work  of  Kitaigoradski .  Equations 
(3)  and  (4)  imply  the  assumed  logarithmic  profile 

U  -  U*/k  log  (Z/Z0)  .  (6) 

Now  the  equation  which  describes  the  magnitude  of 
growth  rate  based  on  the  Miles  instability  theory  and  a  con¬ 
tribution  by  Phillips,  (1966),  can  be  written 


B/f 


pa  2  tt 

(7k 


{ 


A  r2k4 

m 

- — 

cos  0 


U"  -  ,  [U  cos  0  -  C]2  rkZdZ)2 

■  'z* 


+  Ap  0/“  r2  C-U")  cos  0  |U  cos  0  -  C|  £~2kZdZ} 


(7) 


where  p  is  the  density  of  sea  water, 

U"  is  the  mean  wind  profile  curvature, 

IP  is  the  mean  wind  profile  slope, 

0  is  the  directional  difference  between  wind  and  wave, 

k  is  a  constant  and  equal  to  tt, 
m 

Ap  is  a  constant  determined  by  experiment  and  Phillips 
ascertained  that  the  value  is  1.6  *  10  with  an 
uncertainty  of  ±30%, 

T  is  +1  above  the  matched  layer  and  less  than  ♦l 
below  the  matched  layer. 
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The  first  term  on  the  right  hand  side  is  Miles' 
solution  and  the  second  term  on  the  right  side  is  Phillips’ 
contribution.  It  should  be  noted  that  B/f  is  a  dimensionless 
quantity.  By  substituting  the  relationships  in  equations  (3), 
(4),  (5)  and  (6)  in  equation  (7),  it  can  be  shown  that  a  rela¬ 
tionship  exists  between  B/f  and  U*/C  (dimensionless  friction 
velocity).  Several  of  the  above-mentioned  investigators  have 
looked  at  this  relationship  either  on  a  theoretical  basis  or 
in  field  experiments  as  seen  in  Figure  2  (DeLeonibus  and 
Simpson,  1972).  Inoue  investigated  these  results  and  pro¬ 
posed  his  own  curve  which  is  also  shown  in  Figure  2  and  can 
be  defined  as 

B(f,U,)  «  0.39  *  10'V7000  f(U*/C)  -3.1  *  10'2]2 

♦  0.72S  (U*/C)V4  *  10“4(C/U*)2}f  (8) 

If  equations  (2)  and  (8)  were  used  in  a  wave  model 
without  any  constraints,  it  is  conceivable  the  waves  could 
continue  growing  forever.  Obviously,  this  does  not  occur  in 
nature;  therefore,  in  order  to  approximate  the  true  situation, 
the  concept  of  a  fully  developed  sea  was  introduced.  The 
essence  of  this  concept  is  that  if  a  wind  with  the  same  magni¬ 
tude  blows  in  the  same  direction  over  a  given  fetch  for  enough 
time,  the  wave  spectrum  will  become  fully  developed  and,  no 
matter  how  much  longer  the  same  wind  blows,  the  spectrum  will 
no  longer  continue  to  grow.  The  Pierson  and  Moskowitz  fully 
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developed  spectrum  (1964)  was  used  in*,  the  Med  model  and 
has  the  form 

S.M  ■  (9) 

0) 

where  *  2irf 

a  -  8.1  *  10'3 
0  ■  0.74  and 

u>Q  *  g/U^g  5  where  U^g  is  the  wind  speed  as  measured 
19.5  meters  above  sea  level. 

At  each  time  step,  the  spectrum  at  each  grid  point  is  compared 
to  the  Pierson  and  Moskowitz  spectrum  for  the  given  wind.  The 
wave  growth  routine  does  not  allow  the  spectrum  to  exceed  the 
Pierson  and  Moskowitz  spectrum.  This  is  accomplished  by 
modification  of  equation  (1)  into: 

jf  *  [A(S/SJ2}1/2  *  BS](1-(S/SJ2]  ,  (10) 

whose  solution  for  zero  initial  conditions  is: 

S(f,t)  *  A{e.xP  Cgt )  ,  ~  1}  [!  +  {d  exp^ptj  -  1}2j-1/2  .  (u) 

00 

In  order  to  deal  practically  with  the  wave  spectrum,  it  had 
to  be  broken  into  discrete  segments.  At  each  grid  point  the 
total  energy  is  described  by  a  two-dimensional  matrix,  15 
frequency  bands  by  12  direction  bands.  Thus,  for  every  time 
step  180  components  of  the  spectrum  are  investigated  at  each 
grid  point.  Since  there  are  455  sea  points  in  the  Med  model. 
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81,900  pieces  of  information  are  calculated  at  each  time  step. 
The  compass  directions  were  divided  into  30-degree  bands 
because  it  was  assumed  that  the  best  existing  wind  model  could 
not  produce  wind  directions  better  than  ±15  degree  accuracy 
and,  of  course,  the  winds  are  the  driving  force  of  any  wave 
model.  The  frequency  bands  were  limited  to  15  to  speed  the 
computer  calculations  and  reduce  computer  memory  requirements. 
For  14  of  the  frequency  bands,  the  energy  is  summed  over  set 
limits,  but  for  the  fifteenth  band,  the  energy  is  summed 

I 

from  0.164  cycles/sec  to  °°,  which,  of  course,  must  be 
truncated  for  practical  reasons.  The  error  caused  by 
truncation  is  minimal.  Table  1  shows  the  breakdown  of  the 
energy  spectrum  into  the  frequency  and  direction  bands.  The 
fifteenth  frequency  is  always  considered  to  have  a  fully 
developed  spectrum  since  the  highest  frequencies  achieve  the 
fully  developed  spectrum  quite  rapidly. 

Since  equation  (11)  only  provides  a  one-dimensional 
spectrum,  an  equation  developed  by  the  Stereo-Wave  Observation 
Project  (SWOP)  is  used  to  obtain  the  directional  spectra 
(Inoue,  1967).  This  equation  is 

F(w,  6,  U)  =  i  [1  +  (0.50  ♦  0.82  g,”1/2  Cu)U/g)  cos  26 

+  0.32  £'1/2(0)U/g)  cos  4e]  (12) 

for  -ir/2  <  0  <  ir/2  and  0  is  the  angle  between  the  wind  direction 
and  the  wave  direction  and  F(w,  0,  U)  ■  0,  elsewhere. 
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The  directional  spectra  are  computed  for  30°  increments. 
As  an  example,  if  the  wind  direction  were  180°,  the  distribution 
of  the  wave  energy  would  be  as  follows:  37.5%  of  the  energy 

would  be  placed  in  the  180°  direction  band;  25%  would  be 

placed  in  the  150°  and  210°  bands;  and  6.5%  of  the  energy 

would  be  placed  in  the  120°  and  240°  direction  bands.  Of 

course,  if  the  wind  direction  were  190°  rather  than  180°, 
then  the  energy  distribution  would  be  more  biased  in  the  210° 
and  240°  direction  bands  than  in  the  120°  and  150°  bands. 

b.  Wave  Propagation 

Waves  generated  by  the  same  wind  system  tend  to  travel 
as  a  group.  A  wave  train  (a  group  of  waves)  travels  at  a  group 
velocity  which  is  one-half  the  wave  celerity  in  deep  water  and 
exactly  the  wave  celerity  in  shallow  water.  Within  a  group, 
the  various  wave  frequencies  travel  at  different  velocities 
such  that  long-period  (low-frequency)  waves  travel  faster 
than  short-period  (high-frequency)  waves.  Although  waves 
travel  continuously  across  bodies  of  water,  NYU  developed  a 
"jump"  technique  in  the  computer  model  to  propagate  waves. 

With  this  approach,  energy  as  a  function  of  frequency  is 
jumped  from  grid  points  to  grid  points  at  given  time  inter¬ 
vals.  At  every  time  step,  the  energy  growth  is  computed  for 
all  frequencies  at  each  grid  point.  When  the  model  has  deter¬ 
mined  that  enough  time  has  elapsed  for  the  energy  of  a  given 
frequency  to  be  transferred  to  other  grid  points,  then  that 
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energy  is  transferred  all  at  once  instead  of  gradually  as 
would  be  the  case  in  a  gradient  technique.  The  Med  wave 
model  uses  a  conical  conformal  grid,  true  at  40°,  with  a  mesh 
length  of  67  kilometers.  The  grid  system  is  shown  in  Figure  3. 
The  time  step  for  this  grid  system  is  one  hour.  With  these 
grid  constraints,  the  higher  frequencies  are  propagated  every 
four  hours  and  the  lowest  frequencies  are  propagated  every 
hour.  The  propagation  time  for  all  frequencies  is  shown  in 
Table  2.  Although  the  spectral  energies  are  computed  and 
propagated  in  the  NYU  model  on  hourly  increments,  the  "jump" 
technique  does  not  really  depict  a  precise  model  of  the  real 
world  and  some  error  is  expected  in  the  propagation  scheme 
since  energy  does  not  move  in  such  uniform  increments  in 
nature. 

In  the  computer  model,  any  given  grid  point  can  accept 
energy  from  any  of  its  surrounding  grid  points,  depending  on 
the  direction  from  which  the  waves  are  propagated.  Since  the 
propagation  directions  are  divided  into  30°  increments,  some 
difficulty  arises  when  energy  has  to  be  propagated  from  a 
direction  other  than  0°,  90°,  180°  or  270°.  The  problem  is 
depicted  in  Figure  4.  If  the  grid  point  of  interest  is 
depicted  with  an  X  and  the  surrounding  points  are  numbered 
from  one  to  eight,  it  can  be  seen  that  energy  can  be  propagated 
from  points  2,  4,  5  and  7,  which  correspond  to  the  above  direc¬ 
tions  315°,  45°,  225°  and  1358,  respectively.  Thus,  if  the 


waves  are  propagated  W  from  330°,  some  of  the  energy  at  point 
X  will  have  to  come  from  grid  point  1  and  some  energy  from 
grid  point  2.  This  is  accomplished  by  propagating  from  grid 
points  1  and  2  on  alternate  time  steps.  Since  new  wind  fields 
are  read  in  every  six  hours,  both  grid  points  1  and  2  would 
be  examined  on  an  equal  basis. 

c.  Wave  Dissipation 

If  the  wind  is  moving  against  the  waves  while  the 
waves  are  propagating,  the  waves  will  be  dissipated.  Inoue 
developed  a  formula  to  account  for  the  wave  dissipation  if 
the  angle  between  the  wind  and  wave  directions  is  greater 
than  75°.  The  formula  is: 

*  so(f.,0.)[rc^fi4]lc(0i)  (is) 

where  SpCf^.e^  *  spectral  component  after  dissipation, 

S0(fi,6i)  =  spectral  component  before  dissipation 

^i,0i  =  center  frequency  and  direction  of  that  component, 

2 

c  *  constant  690  (for  ft  ), 

S  =  EES  (f . ,8. ) 
w  w v  i  *  1 ' 

and  k(ex)  *  0  for  6  <_  75°, 

k^)  *,1.5  for  0  <  105°, 

k(61)  *  3.0  for  105°  <  0j  <  135°, 

k(01)  -  4.5  for  135°  <  0j  <  165°, 

k(0x)  -  6.0  for  165°  <  e1  <  180°. 
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3.  Mediterranean  Winds 

A  wind  model  which  could  be  used  as  a  driving  force  for 
the  Mediterranean  wave  model  had  to  be  developed.  The  availa 
bility  of  good  wind  analyses  and  prognoses  is  the  key  to 
success  for  any  wave  model.  The  two  wind  models  available 
for  use  on  a  routine  basis  at  FNWC  are  the  Tropical  Band  Grid 
Analysis  with  a  mesh  length  200  kilometers  at  35°N  and  the 
Marine  Wind  Analysis  and  Forecast  model  with  a  mesh  length  of 
400  kilometers  at  35°N.  The  Mediterranean  basin  is  a  rather 
complex  meteorological  area.  Its  many  islands  and  varying 
coastline  create  many  localized  wind  conditions  such  as 
"mistrals"  off  the  coast  of  France  and  "meltemis"  in  the 
Aegean  Sea,  which  cannot  be  detected  in  the  Tropical  Band 
Analysis  or  the  Marine  Wind  Field.  Ideally,  a  detailed  wind 
model  solely  devoted  to  the  Mediterranean  Sea  is  needed.  An 
effort  is  being  made  to  develop  a  new  Med  wind  model;  however 
this  is  a  major  undertaking  and  will  not  be  completed  in  the 
immediate  future.  In  order  to  expedite  the  wave  program,  a 
modified  version  of  the  Marine  Wind  program  was  developed. 

The  Marine  Wind  program  was  chosen  even  though  its  mesh 
length  is  longer  than  the  Tropical  Band  grid  because  it 
produced  both  analyzed  and  forecasted  wind  fields.  The  wind 
field  mesh  length  was  reduced  by  interpolation  from  400 
kilometers  to  67  kilometers.  Ship  reports  were  added  at 
analysis  time  so  that  local  wind  conditions  missed  by  the 
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Marine  Wind  program  would  be  added  to  the  program.  It  was 
recognized  that  the  scheme  would  not  solve  all  the  wind 
definition  problems,  but  it  was  hoped  that  the  errors  would 
be  kept  to  a  minimum. 

It  should  be  noted  that  all  the  wave  growth  formulae 
developed  by  NYU  which  are  a  function  of  wind  speed  or 
frictional  velocity  were  based  on  the  assumption  that  the 
winds  were  19.5  meters  above  the  sea  surface.  The  FNWC 
Marine  Wind  program  also  calculates  the  winds  at  the  19.5 
meter  level;  however,  since  it  would  be  rather  difficult  and 
tedious  to  determine  individual  ships'  anemometer  heights, 
the  wind  velocities  from  the  ship  reports  were  not  corrected 
for  anemometer  height.  Thus,  some  errors  could  be  expected 
in  the  analysis  schemes  since  the  wind  speed  profiles  are 
assumed  to  be  logarithmic  over  the  sea  surface  and  large 
differences  in  wind  velocities  could  exist  over  short  vertical 
distances . 

The  winds  were  added  to  the  wave  model  every  six  hours. 
In  the  North  Atlantic  wave  model,  the  wind  fields  were  added 
at  an  analysis  or  prognosis  time  and  then  used  for  the  next 
six  hours;  that  is,  if  the  wind  input  was  added  at  OZ,  then 
the  same  winds  would  be  used  for  every  time  step  up  to  and 
including  the  05Z  time  step.  It  was  felt  that  the  wind  field 
would  be  more  representative  if  the  winds  were  centered  about 
an  input  time;  that  is,  the  OZ  wind  fields  would  be  used  from 
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21Z  to  02Z  and  the  06Z  fields  would  be  used  from  03Z  to  08Z. 
This  scheme  tended  to  be  superior  to  the  previous  techniques 
except,  as  will  be  discussed  in  Section  4,  when  strong  weather 
systems  entered  an  area  near  the  end  of  an  analysis  or  progno¬ 
sis  period.  Since  the  model  will  use  these  winds  for  the 
entire  six-hour  period,  it  will  produce  large  waves  sooner 
than  they  actually  occur. 
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4.  Mediterranean  Wind  and  Wave  Verification 

The  Mediterranean  wind/wave  model  produces  analyses  and 
forecasts.  Since  it  is  expected  for  obvious  reasons  that 
the  analyses  probably  would  be  consistently  superior'  to  the 
prognoses,  the  verification  section  will  be  divided  into  two 
parts-- (a)  Wind/Wave  Analysis  Verification  and  (b)  Wind/ 

Wave  Prognosis  Verification. 

a.  Wind/Wave  Analysis  Verification 

Although  the  optimum  method  for  verifying  the  wave 
calculations  is  to  compare  the  results  with  data  from  cali¬ 
brated  wave  recorders,  there  is  a  scarcity  of  this  type 
of  data  in  the  Mediterranean.  Thus,  the  model  verification 
was  accomplished  primarily  with  ship  reports.  Data  from  wave 
recorders  will  be  discussed  in  a  later  section. 

(1)  Shipboard  Observations 

Under  most  circumstances  visual  wave  observations, 
especially  by  untrained  observers,  are  not  the  best  source  of 
data;  however,  it  was  hoped  that  the  density  of  ship  reports 
would  be  sufficient  enough  so  that  obvious  incorrect  reports 
could  be  eliminated.  Ship  reports  were  obtained  for  a  time 
period  extending  from  May  to  September  1972.  The  ship  reports 
were  compared  only  when  the  analyzed  winds  were  greater  than 
10  knots.  In  the  original  NYU  North  Atlantic  wave  spectral 
model,  waves  were  not  grown  when  the  wind  speeds  were  less 
than  10  knots.  This  criterion  has  since  been  reduced  to  three 
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knots  so  that  the  energy  in  the  highest  frequency  band  which 
is  not  propagated  can  be  properly  dissipated.  The  difference 
in  wave  heights  generated  by  a  three-knot  wind  and  10-knot 
wind  is  less  than  0.2  meters.  A  total  of  4S1  ship  observa¬ 
tions  was  used.  The  results  for  wind  speed  are  shown  in 
Figures  5  and  6  and  for  wind  direction  in  Figures  7  and  8. 

As  seen  in  Figure  5,  53.7%  of  the  wind  speed  differences 
(analyzed-reported)  were  in  the  ±2.5  knot  category  and  another 
24.8%  fell  into  the  ±5  knot  category.  Less  than  6%  were  in 
categories  greater  than  ±10  knots.  The  reported  wind  speed 
(negative  values)  tended  to  be  greater  than  the  analyzed 
wind  speeds  (Figure  6) ,  particularly  in  the  ±5  knot  range 
where  38.41  of  the  comparisons  were  between  -0.1  knots  and 
-5.0  knots.  In  19.7%  of  the  cases  there  were  no  differences 
at  all. 

The  wind  directions  are  divided  into  30-degree 
bands  (Section  2a).  Over  71.8%  of  the  wind  direction  dif¬ 
ferences  were  in  the  ±15  degree  band  (Figure  7)  and  another 
12%  were  in  the  ±30  degree  band.  Only  8.6%  of  the  differences 
were  greater  than  ±45  degrees.  Figure  8  indicates  that  there 
was  little  error  bias  in  the  wind  direction  comparisons.  It 
would  seem,  based  on  Figures  7  and  8,  that  the  analyzed  wind 
directions  are  quite  good,  but  the  analyzed  wind  speeds  tend 
to  be  slightly  low. 

In  addition  to  shipboard  wind  observations, 
observations  were  also  made  during  airplane  flights  on 
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29-30  May  1972  when  an  airborne  wave  laser  collected  data  for 
model  verification.  These  flights  will  be  discussed  in  greater 
detail  in  a  later  section;  however,  it  is  interesting  to 
compare  the  wind  observations  with  computed  winds  since  the 
airplane  traversed  a  rather  large  area  in  a  relatively  short 
time  and  the  wind  observations  could  be  considered  synoptic. 

The  airplane  flew  at  an  altitude  of  500  feet.  With  knowledge 
of  the  airplane  doppler  navigation  system,  the  measured  wind 
velocities  can  be  reduced  to  ground  level.  The  wind  direc¬ 
tions  should  be  quite  accurate  with  the  measured  speeds 
showing  a  verification  accuracy  to  within  ±5  knots  according 
to  Clinton  F.  Beckner,  NAVOCEANO .  The  airplane  winds  are 
shown  in  Figures  9  and  10.  The  abscissae  of  the  graphs 
indicate  the  time  of  day.  The  geographical  location  of  the 
airplane  at  the  time  can  be  seen  in  Figure  15  .  The  morning 
flight  on  29  May  was  made  over  the  Straits  of  Sicily.  The 
afternoon  flight  on  29  May  and  the  morning  flight  on  30  May 
were  made  southeast  of  Sicily.  The  tail  of  the  arrows  indi¬ 
cates  the  direction  from  which  the  wind  was  coming. 

As  with  the  ship  reports,  the  wind  directions 
compared  more  favorably  than  the  wind  speeds.  Unlike  the  ship 
reports,  the  wind  speed  errors  do  not  seem  to  be  biased  in 
either  a  positive  or  negative  direction. 

If  the  same  criteria  were  used  as  established  for 
the  wind  velocity  comparisons,  then  only  223  wave  height 
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observations  were  available  for  comparison  during  the  same 
time  period.  This  is  not  surprising  since  ships  generally 
do  not  report  waves  as  frequently  as  wind  velocity.  Wave 
directions  cannot  be  compared  as  wind  wave  and  wind  directions 
are  assumed  to  be  the  same. 

The  wave  height  differences  (analyzed-reported) 
are  shown  in  Figures  11  and  12.  Nearly  42.91  of  the  wave 
height  differences  were  within  ±1.5  feet  (approximately  ±0.5 
meters)  and  34.1%  of  the  differences  were  within  ±3  feet  (±1 
meter).  Less  than  9.9%  of  the  differences  were  in  categories 
greater  than  ±4.5  feet.  In  some  of  these  cases,  the  large 
deviations  can  be  explained  by  either  poor  wave  reporting  or 
poor  message  transmission  as  the  ships  reported  wind  speeds 
greater  than  20  knots  but  wave  heights  of  one  foot  or  less. 

The  analyzed  wave  heights  tended  to  be  higher  than  the 
reported  wave  heights.  This  phenomenon  was  not  expected 
since  the  analyzed  wind  speeds  tended  to  be  lower  than  the 
reported  speeds.  Although  a  definite  explanation  of  this 
situation  is  not  available  as  of  this  writing,  several  hypotheses 
exist.  There  is  a  possibility  that  the  wave  energy  in  the  com¬ 
puter  model  grows  at  a  faster  rate  than  it  does  in  the  physical 
world.  The  growth  equations  were  originally  developed  for  the 
North  Atlantic  Ocean  where  the  fetch  and  duration  of  a  given 
wind  are  probably  longer  than  they  would  be  in  the  Mediterranean 
Sea.  It  is  possible  that  the  growth  equations  will  have  to 
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be  refined  for  the  Mediterranean  Sea;  however,  it  is  interesting 
to  note  that  the  only  available  set  of  field  data  where  wave 
heights  increased  from  a  calm  sea  to  approximately  3.0  meters 
indicated  that  the  growth  rate  was  on  the  order  of  0.6  meters 
per  hour  for  a  14  m/sec  wind  speed.  The  computer  model  had 
the  same  growth  rate.  Thus,  the  above  argument  will  have  to 
be  held  in  abeyance  until  additional  data  become  available 
which  can  either  confirm  or  refute  the  argument.  This  data 
set  will  be  discussed  in  more  detail  in  later  paragraphs. 

Another  explanation  may  be  that  the  ship  observers 
underestimate  the  wave  heights.  A  check  was  made  of  the  corre¬ 
lation  between  reported  wind  speeds  and  the  associated  wave 
heights.  It  appeared  that  the  wave  heights  were  much  lower 
than  would  be  expected  for  a  given  wind  speed,  no  matter 
what  the  fetch  might  be.  Generally  speaking,  this  problem 
does  not  seem  to  be  too  serious  since  11%  of  the  wave  height 
differences  fell  between  ±3  feet  (±1  meter). 

Several  decisions  had  to  be  made  before  swell 
heights  and  directions  could  be  compared.  First,  there  has 
been  considerable  speculation  as  to  whether  shipboard 
observers  are  able  to  differentiate  between  wind  waves 
and  swell.  A  check  was  made  of  the  Mediterranean  shipboard 
observations  to  determine  if  observers  were  dif ferentinating 
between  wind  waves  and  swell.  It  was  determined  that  in  641 
of  the  ship  reports  evaluated,  the  swell  and  wind  (wave 
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direction  assumed  to  be  the  same  as  wind)  were  within  +60 
degrees  of  each  other.  If  the  climatological  and  geographical 
make-up  of  the  Mediterranean  basins  is  considered,  this  ratio 
probably  is  not  out  of  proportion;  however,  it  was  decided 
to  evaluate  the  reported  swell  directions  very  carefully 
before  comparing  the  data  to  the  model  results. 

Another  point  to  consider  is  the  method  by  which 
the  Mediterranean  wave  model  computes  wave  heights  and  direc¬ 
tions.  As  mentioned  in  Section  2a,  the  model  computes  and 
propagates  wave  energy  rather  than  wave  heights.  At  any 
given  grid  point,  the  wave  energy  is  the  sum  total  of  all  the 
energy  in  the  12  direction  bands  by  15  frequency  bands  matrix. 

It  would  be  impractical  in  terms  of  computer  output  time  and 
the  volume  of  output  generated  to  display  the  spectra  at 
every  grid  point.  Instead,  the  spectra  are  displayed  at 
representative  grid  points  and  significant  wave  heights, 

^1/3  (ky  definition  the  average  height  of  the  highest  one- 
third  waves)  is  displayed  for  all  grid  points.  The  signifi¬ 
cant  wave  height  by  definition  is 

**1/3  *  4/ET0TAL  (14) 

where  E-roTAL  is  the  total  wave  energy  at  a  grid  point.  The 
wave  direction  is  determined  by  investigating  the  total  energy 
in  each  direction  band  and  then  selecting  the  band  which  con¬ 
tains  the  modal  energy.  If  the  modal  wave  direction  at  a  grid 
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point  is  not  the  same  as  the  wind  direction,  then  the  wave 
is  assumed  to  be  swell  rather  than  a  wind  wave.  One  problem 
that  arises  with  this  approach  is  that  if  a  secondary  wave 
front  exists,  it  could  be  difficult  to  ascertain  from  a 
printout  of  wave  heights  and  directions.  The  spectra  output 
for  the  selected  grid  points,  which  hopefully  represent  the 
entire  Mediterranean  Sea,  would  have  to  be  carefully  analyzed. 

Because  of  the  above  limitations,  it  was  decided 
to  use  only  ship  data  from  10-15  April,  24-28  April  and 
25-30  May  1972  for  swell  comparisons.  These  dates  were 
selected  because  hindcast  studies  had  been  made  with  the  wave 
spectral  model  for  comparison  with  data  from  a  SACLANT 
(Supreme  Allied  Command- Atlantic)  ASW  Research  Centre  wave 
buoy  and  NAVOCEANO  airborne  wave  laser.  The  analyses  and 
ship  reports  for  these  time  periods  were  reviewed  very  care¬ 
fully.  Comparisons  were  only  made  when  it  was  felt  that  the 
ships  were  actually  reporting  swell  and  not  wind  waves. 

Even  with  careful  quality  control,  there  were  still  ship 
reports  which  were  used  for  comparison  that  might  have  been 
eliminated. 

The  comparisons  between  analyzed  and  reported 
swell  heights  and  directions  are  plotted  in  Figures  13  and  14, 
respectively.  Since  approximately  30%  of  the  ship  reports  con¬ 
tained  swell  observations  and  many  of  these  reports  were  not 
available  because  of  the  above  restraints,  only  61  ship  reports 
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were  available  for  comparison.  Over  901  of  the  differences 
(analyzed-reported)  were  within  ±3  feet  (±1  m) ,  with  70.51 
within  ±1.5  feet  (±0.5  m) .  The  swell  directions  compared 
almost  as  well  as  the  swell  heights.  About  72.1%  of  the 
differences  were  within  ±30  degrees  and  another  14.7%  within 
±60  degrees.  Neither  the  swell  height  or  direction  dif¬ 
ferences  seemed  to  be  biased  in  either  a  positive  or 
negative  direction. 

Some  of  the  direction  differences  could  be  readily 
explained.  In  one  case  where  the  error  was  +180  degrees,  it 
was  quite  apparent  that  the  computer  model  took  approximately 
three  hours  longer  than  it  should  have  to  propagate  swell 
which  was  coming  from  due  west  to  the  ship  position  (33.9°N, 
28.3°E).  Since  there  was  only  one  ship  report  from  all  three 
studies  available  in  the  far  eastern  Mediterranean  Sea,  it  is 
difficult  to  determine  if  the  model  has  an  inherent  propagation 
problem  in  the  eastern  Mediterranean  Sea,  where  the  swell 
can  move  in  a  longer  continuous  path  than  in  any  other  area 
in  the  Mediterranea,  or  if  this  was  an  isolated  problem.  This 
situation  does  not  seem  to  occur  in  the  western  or  central 
Mediterranean  basins;  however,  in  these  areas,  swell  can  only 
propagate  over  relatively  short  distances  before  encountering 
land.  In  at  least  two  other  cases  where  the  differences  were 
equal  to  or  greater  than  ±90  degrees,  computed  swell  probably 
did  exist  at  the  ship  position  but  the  printout  of  wave  spectra 
from  the  closest  grid  points  seemed  to  indicate  that  the  computed 
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wind  waves  dominated  the  computed  swell  rather  than  vice  versa 
as  indicated  by  the  ship  reports.  Clearly  this  type  of  problem 
can  only  be  resolved  with  wave  data  measured  by  accurate  recorders. 
Some  of  the  differences  (of  between  +30  and  60  degrees)  can  be 
attributed  to  the  method  by  which  the  computer  model  propagates 
waves.  As  discussed  in  Section  2b,  the  computer  model  does  not 
propagate  wave  energy  from  all  compass  directions  at  every 
time  step;  instead  it  propagates  from  alternate  directions. 

Thus,  at  any  given  time  step,  the  computed  wave  direction  could 
vary  by  as  much  as  30  degrees  from  the  previous  time  step.  If 
one  considers  all  the  assumptions  that  had  to  be  made  in  invesgi- 
gating  the  swell  parameters  and  that  the  Mediterranean  Sea 
contains  many  little  islands  which  could  easily  obstruct  swell 
propagation,  then  it  does  appear  that  the  model  handles  swell 
propagation  satisfactorily.  Waves  will  again  be  analyzed  in 
the  next  section  in  terms  of  frequency  or  period.  Since 
measured  wave  data  will  be  reviewed,  fewer  assumptions  will 
have  to  be  made  than  were  made  with  shipboard  observations. 

(2)  Measured  Wave  Data 

The  only  measured  wave  data  available  for  this 
study  were  obtained  from  a  Datawell  Waverider  Buoy  System 
operated  by  the  SACLANT  ASW  Research  Centre  and  an  airborne 
wave  laser  operated  by  NAVOCEANO.  Measurements  from  the 
Waverider  were  made  on  10-15  April  and  25-29  May  1972.  On 
the  first  occasion  the  buoy  was  located  at  43°30'N,  08°30fE 
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(Ligurian  Sea).  These  data  were  only  compared  to  the  computer 
wave  analysis.  On  the  second  occasion  the  buoy  was  located 
at  37°23'N,  11°36'E  (Straits  of  Sicily).  These  data  were  com¬ 
pared  to  both  computer  analyses  and  24-hour  prognoses.  The  buoy 
locations  are  shown  in  Figure  15.  The  Datawell  Waverider 
actually  measures  vertical  acceleration  and  must  be  doubly 
integrated  to  obtain  wave  heights.  The  Datawell  Waverider 
Buoy  System  has  proved  to  be  a  highly  reliable  and  accurate 
instrument  for  measuring  waves.  The  wave  spectra  computed 
from  the  Waverider  were  supplied  by  Dr.  Melbourne  Briscoe, 
SACLANT  oceanographer. 

The  computed  and  measured  wave  spectra  cannot  be 
compared  exactly  since  the  wave  spectra  from  the  buoy  data  are 
computed  for  individual  frequencies  while  the  wave  energies  in 
the  computer  model  are  computed  for  frequency  bandwidths.  The 
distribution  of  the  bandwidths  is  shown  in  Table  2.  First, 
it  should  be  noted  that  the  bandwidths  are  not  equal  in  size. 
This  tends  to  produce  some  unevenness  in  the  wave  spectra. 

The  high  frequency  band  theoretically  includes  all  the  gravity 
wave  energy  from  0.164  to  00  cycles/sec.  Thus,  the  first  point 
in  the  computed  wave  spectra  usually  appears  as  a  spike.  To 
compare  the  computed  and  measured  wave  spectra  exactly,  the 
wave  energy  for  the  individual  frequencies  in  the  measured 
spectra  would  have  to  be  summed  and  then  compared  to  the 
energy  in  a  given  frequency  band  of  the  computed  spectra. 
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Since  this  would  be  a  rather  arduous  task,  two  other  parameters 
that  can  be  more  easily  deduced  from  the  spectra  will  be 
compared.  These  parameters  are  the  dominant  frequencies  and 

fil/3‘ 

The  measured  wave  spectra  from  13-14  April  were 
compared  to  computed  data  from  grid  points  15  (43.89°N, 

7.96°E)  and  26  (43.23°N,  8.01°E)  and  are  shown  in  Figures 
16-21.  Since  there  was  an  interval  of  calm  seas  during 
25-29  May,  only  selected  wave  spectra  are  shown  in  Figures 
24-29.  These  spectra  demonstrate  the  growth  of  waves  from 
calm  to  turbulent  conditions.  During  this  time  period,  com¬ 
puted  spectra  from  grid  points  220  (37.32°N,  11.82°E),  221 
(37 . 33°N,  12 . 66°E) ,  250  (36.55°N,  11.83°E)  and  251  (36.66°N, 
12.67°E)  were  used  for  comparison.  The  laser  wave  spectra 
comparisons  are  shown  in  Figures  30-31.  Generally  speaking, 
the  computed  spectra  compared  rather  well  to  the  measured 
spectra.  Some  variation  did  occur  in  several  of  the  compari¬ 
sons  and  these  anomalies  v/ill  now  be  described. 

First,  the  significant  wave  height  differences  for 
the  13  April  wave  spectra  (Figures  16-19)  were  less  than  0.1 
meters.  On  14  April  (Figures  20-21)  the  differences  increased 
to  approximately  0.5  meters.  The  source  of  error  probably  can 
be  attributed  to  the  14  April,  12Z  computer  wind  analysis  in 
the  vicinity  of  the  wave  buoy  implantment.  The  wind  analysis 
was  strongly  influenced  by  three  ship  reports  that  contained 
wind  speeds  of  25  knots  and  were  120  to  170  nautical  miles 
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(three  or  four  grid  mesh  lengths)  west  of  the  buoy.  The  true 
wind  speeds  at  the  buoy  site  probably  were  less  than  the  cal¬ 
culated  wind  speeds;  however,  the  SACLANT  ASW  Research  Centre 
did  not  provide  wind  data  for  this  time  period  and  there  were 
no  other  ships  in  the  immediate  area  to  provide  the  actual 
wind  speeds  in  the  area  of  interest.  It  is  interesting  to 
note  that  the  second  wave  spectrum  based  on  14  April  buoy 
data  did  increase  slightly  in  magnitude,  but  since  the  time 
sequence  did  not  extend  any  further,  it  could  not  be  determined 
if  the  measured  H^^  would  eventually  approach  the  computed 
^1/3*  A  similar  situation  occurred  on  28  May  when  the  computed 
Hj/j  at  12Z  was  almost  twice  as  large  as  the  measured  H^y^ 
(Figures  25-26)  and  two  hours  later  the  difference  in  magnitude 
of  the  significant  wave  heights  was  approximately  0.1  meter. 
This  situation  is  definitely  a  function  of  the  method  used  by 
the  model  to  input  winds  (Section  3).  This  set  of  data  will 
be  discussed  further  in  the  next  several  paragraphs. 

On  13  and  14  April,  the  dominant  frequencies  of 
the  measured  data  were  on  the  order  of  0.2  Hz.  The  dominant 
frequency  band  of  the  computed  spectra  was  the  highest  band 

(0.164  Hz  to  °°) .  Thus,  there  is  a  high  degree  of  correlation 

% 

between  the  computed  and  measured  dominant  frequencies. 

Dr.  Briscoe  states  that  the  nature  of  waves  in 
the  western  Mediterranean  Sea  is  such  that  the  predominant 
frequency  usually  lies  between  0.2  Hz  and  0.3  Hz.  The  computer 
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model  does  not  presently  have  this  degree  of  resolution,  but 
the  computed  high  frequency  band  could  be  shifted  to  account 
for  the  higher  frequencies  at  the  expense  of  the  lower  fre¬ 
quency  bands  which  apparently  have  little  significance  in  the 
Mediterranean  Sea.  An  unpublished  climatology  study  (Lazanoff, 
1972)  indicates  that  the  percentage  of  wave  frequencies  less 
than  0.05  Hz  was  on  the  order  of  2%  or  less. 

The  first  two  computed  spectra  on  13  April  had  a 
second  energy  peak  in  the  vicinity  of  0.1  Hz.  The  measured 
data  did  not  have  this  peak.  Since  the  secondary  energy  peak 
significantly  diminished  by  18Z,  13  April,  this  may  be  an 
indication  that  the  wave  model  is  not  dissipating  energy  as 
fast  as  it  should.  Again,  since  there  was  only  one  set  of 
data  where  this  situation  occurred,  it  is  difficult  to  make 
a  general  statement  about  the  dissipation  techniques  of  the 
computer  model. 

For  the  most  part,  the  wave  spectra  comparisons 
for  25-29  May  seem  to  be  better  than  the  April  comparisons. 

Wind  speeds  and  H-^  are  plotted  for  the  entire  time  period 
in  Figures  22  and  23.  As  mentioned  previously,  the  computed 
Hi/3  was  almost  twice  as  large  as  the  measured  at  12Z, 

28  May  (Figure  25)  and  two  hours  later,  the  difference  in 
significant  wave  heights  was  on  the  order  of  0.1  meter 
(Figure  26).  The  rather  large  error  at  12Z  can  be  attributed 
to  the  wind  cycling  method  used  by  the  model  to  input  wind 
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velocities.  This  can  be  demonstrated  by  examining  the  graphs 
in  Figure  23.  The  computed  and  measured  wave  heights  and  wind 
speeds  from  20Z,  27  May  to  12Z,  29  May  are  plotted.  Wind 
directions  were  not  shown  since  the  directions  never  varied 
from  the  300°  to  330°  range  during  the  entire  time  period. 
Although  the  plots  indicate  there  is  good  agreement  between 
the  wind  analyses  and  measured  wind  speeds,  it  must  be 
remembered  that  the  wind  analyses  are  only  produced  every 
six  hours  and  cannot  reflect  any  fast-moving  changes  such  as 
occurred  between  04-06Z,  28  May.  Thus,  if  the  winds  have 
high  velocities  for  a  short  duration  near  analysis  time,  the 
model  could  over-grow  waves  since  the  same  wind  analysis  is  used 
for  six  hours  CSection  3) .  This  is  precisely  what  occurred  at 
12Z,  28  May.  Although  the  wind  velocities  did  not  become  large 
until  11Z,  the  wind  analysis  reflected  these  winds  at  09Z, 
permitting  the  model  to  grow  the  large  waves  two  hours  sooner 
than  they  actually  occurred.  The  inverse  could  also  occur, 
i.e.,  the  analyzed  winds  could  be  too  low  for  several  hours 
and  the  wave  growth  would  lag  by  the  same  amount  of  time.  The 
obvious  solution  to  this  problem  would  be  to  increase  the 
number  of  input  wind  fields  per  time  steps.  FNWC  is  investi¬ 
gating  a  more  sophisticated  wind  analysis  for  the  Mediterranean. 
The  wind  fields  will  be  produced  on  a  fine  mesh  length  and  be 
based  on  both  land  and  sea  reports.  If  the  product  is  better 
than  the  present  wind  analysis,  then  the  wind  analysis  tech¬ 
niques  used  in  the  wave  spectral  model  will  be  modified. 
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As  of  now  the  analysis  time  can  be  considered  accurate 
within  ±3  hours. 

The  wave  heights  shown  in  Figure  23  indicate  that, 
given  the  correct  wind  speed,  the  model  appears  to  be  growing 
waves  at  the  proper  rate.  This  subject  has  already  been  dis¬ 
cussed  in  the  paragraphs  describing  the  analysis  of  observed 
wave  heights  from  ships. 

The  airborne  wave  laser  was  installed  aboard  a 
NAVOCEANO  P-3A  airplane.  Flights  were  made  on  29-30  May  1972. 
The  wave  measuring  instrument  is  a  standard  Geodolite  3-A  air¬ 
borne  altimeter,  which  is  a  continuous  wave  helium-neon  laser, 
manufactured  by  Spectra  Physics,  Inc.  The  helium-neon  laser 
is  described  in  detail  by  Ross,  Peloquin  and  Sheil  (1968)  and 
Ross,  Cardone  and  Conaway  (1970).  The  airplane  tracks  were 
flown  in  pairs--upwind  and  downwind- -at  an  altitude  of  500 
feet.  This  approach  is  used  so  that  the  speed  of  the  aircraft 
relative  to  the  phase  speed  of  each  wave  frequency  component 
can  be  accounted  for  without  too  much  difficulty  when  the  moving 
coordinate  system  of  the  aircraft  is  converted  to  fixed  coordi¬ 
nates.  Since,  in  this  case,  it  is  assumed  that  all  waves  are 
traveling  in  the  direction  of  the  wind,  the  presence  of  swell 
from  other  directions  and  the  spreading  of  wave  energy  with 
distance  can  lead  to  errors  [Schule,  Simpson  and  DeLeonibus 
(1971)].  During  25-29  May  the  range  of  wind  directions  remained 
within  280°-320°  so  that  the  swell  was  moving  in  the  same 
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direction  as  the  wind  waves  and  this  type  of  error  should  be 
at  a  minimum.  Errors  due  to  wave  dispersion  should  also  be 
at  a  minimum  since  airplane  tracks  for  29  May,  as  shown  in 
Figure  15,  were  made  in  narrow  areas.  Unfortunately,  due  to 
instrument  failure,  the  laser  data  on  30  May  were  inaccurate 
and  will  not  be  discussed  in  this  paper. 

Airplane  motion  can  distort  the  laser  data.  The 
full  effect  of  the  motion  is  not  fully  understood.  Ross  et 
al.  (1970)  indicated  that  the  disturbances  would  probably 
occur  in  the  low  frequency  range  (.065  Hz  or  less).  Since  a 
wave  climatology  study  (Lazanoff,  unpublished  report)  has  indi¬ 
cated  that  the  swell  frequencies  seldom  reach  or  exceed  this 
range,  errors  due  to  airplane  motion  should  also  be  at  a 
minimum. 

Usually  waves  are  measured  by  a  stationary  wave 
recorder  and  wave  spectra  are  computed  from  a  finite  length 
of  data.  Since  the  airplane  is  not  stationary,  the  computa¬ 
tion  of  wave  spectra  from  the  laser  is  complicated.  The  air¬ 
plane  flew  at  speeds  that  ranged  from  235  to  265  knots.  Thus, 
the  laser  wave  spectra  is  a  function  of  distance  as  well  as 
time.  If  a  small  enough  time  interval  can  be  selected  and  the 
airplane  flies  over  deep  water,  then  it  can  be  assumed  that  the 
wave  spectrum  is  representative  of  one  location.  In  this  case, 
the  laser  recorded  data  every  second  for  30  minutes  on  each 
track.  The  spectra  were  computed  for  three-minute  intervals 
(approximately  12,5  miles). 
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Comparisons  of  spectra  from  the  wave  laser  with 
the  computed  spectra  are  shown  in  Figures  30  and  31.  Spectrum 
from  the  wave  buoy  is  also  shown  in  Figure  30.  The  wave  buoy 
spectrum  agrees  with  the  laser  spectrum.  The  significant 
wave  heights  and  the  predominant  frequencies  (approximately 
0.13  Hz)  compared  extremely  well.  Both  measured  spectra 
show  a  secondary  frequency  peak  at  0.1  Hz  which  does  not 
appear  in  the  computed  spectrum.  The  secondary  energy  peak 
is  a  definite  indication  of  swell.  Since  the  model  does  not 
have  this  peak,  this  may  be  an  error  similar  to  the  problem 
noted  in  the  swell  analysis  section  where  it  seemed  that  the 
swell  did  not  propagate  energy  fast  enough  in  the  eastern 
Mediterranean. 

The  comparison  between  ship  reports,  measured 
data  and  computer  calculations  seems  to  indicate  that  the 
wind  and  wave  analyses  have  a  high  degree  of  accuracy.  The 
wave  analyses,  as  well  as  the  24-hour  prognoses,  are  being 
transmitted  to  Rota,  Spain,  twice  a  day  (0Z  and  12Z)  on  an 
operational  basis. 

b.  Wind/Wave  Prognosis  Verification 

Since  it  has  been  shown  in  the  previous  section  that, 
given  the  correct  wind  velocity,  the  computer  model  will  pro¬ 
duce  accurate  wave  spectra,  the  principal  concern  of  this 
section  is  to  determine  if  the  wind  prognoses  are  valid.  To 
reiterate  what  has  already  been  mentioned  in  Section  3: 


the  Med  wind  prognoses  are  interpolated  from  the  FNWC  Marine 
Wind  fields.  The  grid  of  the  Marine  Wind  model  has  a  mesh 
length  on  the  order  of  400  kilometers  at  35°N.  The  Mediterranean 
grid  has  a  mesh  of  67  kilometers;  thus,  most  local  wind  condi¬ 
tions  will  not  be  predicted  with  this  scheme.  Unlike  the  Med 
wind  analyses,  the  prognoses  cannot  be  modified  with  ship 
reports.  Wind  prognoses  are  made  every  six  hours  out  to 
48  hours.  It  can  be  assumed  that  the  more  time  that  has 
elapsed  from  the  analysis  time,  the  less  accurate  the  prog¬ 
nosis.  Since  FNWC  decided  to  transmit  the  24  prognoses  to 
FWC,  Rota,  only  these  prognoses  are  discussed  in  this  section. 

Wind  prognoses  from  three  grid  points--175  (37.5°N, 

4 . 5°E) ,  220  (37 . 5°N,  12°E),  and  372  (3S°N,  28°E)--were 
selected  for  verification.  The  grid  points  were  to  repre¬ 
sent  the  Eastern  and  Western  Mediterranean  Basins  and  the 
Straits  of  Sicily.  The  prognoses  for  grid  points  175  and 
220  were  compared  to  ship  reports  when  available  (there  were 
no  ship  reports  available  in  the  vicinity  of  point  372  during 
the  verification  period) .  Prognoses  for  grid  point  220  were 
also  compared  to  measured  winds  in  the  vicinity  of  the  SACLANT 
wave  buoy  for  the  period  25-29  May.  For  the  remainder  of  the 
verification  period  (30  May-6  June),  the  prognoses  were  com¬ 
pared  to  the  computer  analyses  and  the  FWC,  Rota  24-hour 
prognoses  (26  May-1  June).  The  same  comparisons  were  made 
during  28  May-6  June  for  grid  points  175  and  372.  In  addition, 
24-hour  prognoses  were  compared  to  wind  velocities  measured  by 
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the  NAVOCEANO  airplane  on  12Z,  29  May  and  12Z,  30  May.  The 
above  comparisons  are  shown  in  Figures  32-36.  The  plotted 
wind  directions  are  the  directions  toward  which  the  winds 
are  moving. 

Although  there  were  not  many  cases  of  high  wind 
speeds,  one  significant  observation  can  be  made  from  the  plots. 
The  prognoses,  unlike  the  analyses,  tend  to  lag  behind  the 
actual  movement  of  weather  systems  by  approximately  12  hours. 
This  was  not  totally  unexpected.  The  Marine  winds  are  computed 
from  surface  pressure  fields  generated  by  the  FNWC  Primitive 
Equation  (PE)  multi-layer  model.  It  has  been  shown  that  the 
PE  model  has  a  tendency  to  move  weather  systems  too  slowly 
over  the  North  Pacific  by  15%  and  over  the  North  American 
continent  by  25%  (Osburn,  1971).  It  would  seem  reasonable 
that  this  tendency  would  extend  to  the  Mediterranean  Sea. 

The  computer  prognoses  were  compared  to  the  FWC, 

Rota  hand  prognoses  to  determine  if  one  method  was  consis¬ 
tently  superior  to  the  other.  The  computer  prognoses  appeared 
to  be  superior  at  grid  points  220  and  372.  For  example,  as 
seen  in  Figure  1,  the  hand  prognosis  completely  missed  the 
weather  system  which  passed  over  grid  point  220  at  approxi¬ 
mately  12Z,  28  May.  The  computer  prognosis  did  predict  the 
system,  although  12  hours  too  late.  The  hand  prognosis 
seemed  to  do  better  at  grid  point  175.  The  superiority  may 
be  due  to  the  geographical  proximity  of  FWC,  Rota  to  the 
grid  point.  Several  U.S.  Navy  air  squadrons  are  stationed 


at  Rota  and  the  weather  central  may  be  able  to  obtain 
meteorological  information  from  these  squadrons  which  is 
not  available  to  anyone  else.  Although  there  wasn’t  much 
difference,  the  computer  prognoses  did  seem  to  compare 
better  to  the  airplane  winds  as  seen  in  Figure  36  than  the 
hand  prognoses. 

Another  interesting  point  about  the  computer 
prognoses  that  was  not  revealed  in  the  above  study  is  that 
in  the  northwest  Mediterranean  (in  the  vicinity  of  France), 
a  consistent  error  seems  to  occur.  The  model  tends  to  over¬ 
compute  wind  speeds  in  this  area.  This  may  be  due  to  the 
close  proximity  of  the  French  Alps  to  the  Mediterranean. 

Although  pressure  fields  over  mountains  are  supposed  to  be 
corrected  to  sea  level  pressure,  the  predicted  pressure 
gradients  seemed  to  be  rather  tight.  This  may  be  due  to 
the  coarseness  of  the  original  grid  or  the  inaccuracy  of  the 
sea  level  pressure  reduction. 

Although  the  Med  wave  spectral  model  has  been 
operational  since  April  1972,  significant  storms  did  not 
occur  in  the  Mediterranean  until  late  October  1972.  Twenty- 
four-hour  forecasts  and  analyses  of  wave  heights  are  shown 
in  Figure  37  for  31  October- 2  November.  Ship  reports 
are  indicated  on  the  analysis  plots.  The  analyses,  prog¬ 
noses  and  ship  reports  show  good  agreement.  During  this  time 
period,  a  storm  passed  over  the  Aegean  Sea.  Maximum  significant 
wave  heights  were  on  the  order  of  15  feet.  The  analyzed  and 
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forecasted  wind  velocities  are  plotted  in  Figure  38  for  a 
selected  grid  point,  194,  in  the  Aegean  Sea.  The  forecasts 
and  the  analyses  compared  quite  satisfactorily  with  each 
other  as  well  as  with  ship  reports.  As  of  12Z,  31  October, 
the  prognoses  appeared  to  lag  behind  the  analysis  by  12  hours. 
The  lag  can  be  seen  more  easily  in  the  wind  direction  plots 
than  in  the  wind  speed  plots.  The  movement  of  the  wave 
system  shown  in  Figure  37  seems  to  confirm  the  prognosis  lag. 
It  may  be  that  the  rate  of  movement  of  forecasted  weather 
systems  in  the  Mediterranean  will  have  to  be  adjusted  by 
artificial  means. 
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5.  Conclusions 

The  Mediterranean  wave  spectral  computer  model  is  the 
first  spectral  model  to  be  used  operationally  for  analyses 
and  prognoses.  The  model  has  two  basic  parts- -wind  computa¬ 
tion  and  wave  energy  computation.  The  wind  analyses  were 
compared  to  ship  reports  and,  for  a  two-day  period,  winds 
were  measured  from  an  airplane.  Although  the  analyzed  wind 
speeds  ind  wind  directions  compared  rather  well  with  the 
observed  data,  the  wind  directions  appeared  to  be  more  accu¬ 
rate  than  the  wind  speeds  which  tended  to  be  too  low.  For 
the  most  part,  the  analysis  scheme  seemed  to  propagate 
weather  systems  at  the  proper  rate  of  speed;  however,  some 
error  did  seem  to  occur  because  the  wind  field  inputs  to  the 
wave  model  were  only  changed  every  six  hours  Cone  hour  per 
time  step).  Although  further  investigation  is  needed  on  this 
subject,  it  would  seem  that  if  the  wind  inputs  were  increased 
to  one  every  three  time  steps,  significant  improvement  in  the 
wind  analysis  would  occur.  FNWC  is  currently  investigating 
a  new  wind  analysis  scheme  for  the  Mediterranean  Sea.  If 
the  new  scheme  proves  to  be  a  significant  improvement  over 
the  current  operational  wind  analysis,  then  it  will  be 
inserted  into  the  Mediterranean  model. 

The  wind  prognoses  were  not  compared  to  measured  data 
as  extensively  as  the  wind  analysis;  however,  it  was  apparent 
from  observation  of  the  24-hour  prognoses  that  the  prognoses 
generally  lagged  behind  the  analysis  by  12  hours.  This 
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error  could  be  attributed  to  the  surface  pressure  fields 
computed  by  the  FNWC  PE  boundary  layer  model  which  are  the 
initial  input  to  the  prognosis  scheme.  The  wind  prognoses 
could  also  be  improved  by  developing  a  finer  mesh  length  for 
the  initial  wind  calculations . 

The  computed  wave  spectra  and  significant  wave  heights 
were  compared  to  ship  reports,  wave  spectra  obtained  from  a 
SACLANT  Wave  Buoy  and  a  NAVOCEANO  airborne  wave  laser.  The 
basic  conclusion  is  that,  given  the  correct  wind  velocity, 
the  computer  model  computes  the  proper  spectra.  It  is 
interesting  to  note  that  the  wave  energy  growth  equations 
which  were  originally  developed  for  the  North  Atlantic  Ocean 
seemed  to  be  accurate  in  the  Mediterranean  Sea.  Two  points 
could  be  investigated  in  greater  detail.  First,  the 
frequency  bandwidths  of  the  computed  wave  spectra  could  be 
modified  so  that  the  highest  frequency  band  is  subdivided 
into  smaller  bandwidths.  This  would  not  affect  the  basic 
calculations  of  wave  energy,  but  would  permit  closer  obser¬ 
vation  of  the  energy  levels  in  the  higher  frequency  ranges. 
This  could  be  of  use  to  underwater  acousticians  who  are 
interested  in  high-frequency  ambient  noise  problems. 

Finally,  a  limited  amount  of  data  seemed  to  indicate  that 
the  model  was  propagating  swell  (or  energy  in  the  lower 
frequencies)  too  slowly.  A  field  test  should  be  conducted 
to  determine  if  the  swells  are  being  propagated  at  the  proper 
rate . 
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Central  Central  Frequency 

Frequency  (H  )  Period  (Seconds)  Bandwidth  (H  ) 

Z  z 


0.164 

6.1 

.164 

- 

00 

0.153 

6.5 

.142 

- 

.164 

0.133 

7.5 

.125 

- 

.142 

0.117 

8.6 

.108 

- 

.125 

0.103 

9.7 

.097 

- 

.108 

0.092 

10.9 

.086 

- 

.097 

0.083 

12.0 

.080 

- 

.086 

0.078 

12.9 

.075 

- 

.080 

0.072 

13.8 

.069 

- 

.075 

0.067 

15.0 

.064 

- 

.069 

0.061 

16.4 

.058 

- 

.064 

0.056 

18.0 

.053 

- 

.058 

0.050 

20.0 

.047 

- 

.053 

0.044 

22.5 

.042 

- 

.047 

0.039 

25.7 

.036 

- 

.042 

Energy  in  these  frequency  bands  is  ordered  in  twelve  (12) 
30  degree  direction  bands  starting  with  the  cardinal  direction 
north 


TABLE  1 
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MEDITERRANEAN  SPECTRAL  FREQUENCIES  VS  DISTANCE  TRAVELED  /HR. 
AND  TOTAL  TIME  FROM  GRID  POINT  TO  GRID  POINT 


DISTANCE  (NM)  ELAPSED  TIME  (HR) 

FREQUENCY  (  HZ)  TRAVELED  /HOUR  FROM  GRID  POINT  TO 

GRID  POINT 


0.15278 

9.75 

4 

0.13333 

II  .25 

3 

0.11667 

12.90 

3 

0.10278 

14.55 

2 

0.09167 

16.35 

2 

0.08333 

18.00 

2 

0.07778 

19  35 

2 

0.072222 

20.7 

1 

0.06667 

22  5 

1 

0.061  II 

24.6 

1 

0.05556 

27.0 

1 

0.05000 

30.0 

1 

0.04444 

33.75 

1 

0.03889 

38.55 

1 

TABLE  2 
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Wind  Velocity  ( knot*  of  6.1m ) 


FIGURE  1 


NUMBER  OF  OBSERVATIONS 
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WINO  SPEED  DIFFERENCES 
FOR  WIND  -  SPEEDS  GREATER 
THAN  10  KNOTS 
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WIND  DIRECTION  DIFFERENCES 
FOR  WIND  SPEEPS  GREATER 


FIGURE  8 


29  MAY  1972 


TIME  (GMT) 

FIGURE  9 
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WAVE  HEIGHT  DIFFERENCES  IN  FEET 


WAVE  HEIGHT  DIFFERENCES  FOR  WIND 
SPEEDS  GREATER  THAN  10  KNOTS 
ANALYZED  WAVE  HEIGHTS  -  REPORTED 


FIGURE  12 
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SACLANT  WAVE  BOUY  EMPLANTMENT  SITES  AND  NAVOCEANO  AIRPANE  TRACKS 
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ANALYZED  WAVE  SPECTRA 

-  COMPUTER  MODEL 
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1800  Z 

H  1/3  *  1.22  M  (  4.00  FT) 

-  SACLANT  WAVE  DATA 
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ANALYZED  WAVE  SPECTRA 


-  COMPUTER  MODEL 
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ANALYZED  WAVE  SPECTRA 
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ANALYZED  WAVE  SPECTRA 
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COMPARISON  OF  COMPUTED  AND  MEASURED  WIND  VELOCITIES 
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COMPARISON  OP  WIND  ANALYSES  AND  PROGNOSES 
AT  37°  30'N,  I2°E 
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COMPARISON  OF  WIND  ANALYSES  AND  PROGNOSES 
AT  37.5°  N,  4.5°  E 
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COMPARISON  BETWEEN  MEASURED  WINDS  FROM  AIRPLANE 
3  15-4-  AND  WIND  PROGNOSES 
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COMPARISON  OF  WIND  ANALYSES  AND  PROGNOSES  (37.32°N,  26.I5°E) 
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The  Mediterranean  wave  model  uses  a  conic  conformal  grid,  permitting  the 
assumption  of  equal  spacing  between  grid  points.  There  are  455  sea 
points  with  a  mesh  length  of  67  km. 

The  wave  model  driving  force  is  a  modified  version  of  the  FNWC  Marine 
Wind  model  which  has  a  mesh  length  of  approximately  370  kms .  For  use  in 
the  Mediterranean,  the  winds  are  interpolated  between  Marine  Wind  field 
grid  points.  At  analysis/pre- analysis  times,  wind  reports  from  synoptic 
ship  files  are  reanalyzed  in  the  Mediterranean  wind  program  so  that  local 
wind  phenomena,  such  as  ^mistrals,*  are  included  in  the  wind  field.  r 

The  Mediterranean  wind  and  wave  model  was  evaluated  with  ship  reports, 
one  wave  buoy  and  an  airborne  wave  laser.  The  wave  computations  appear 
to  be  as  good  as  the  wind  input  allows. 

The  on-time  analyses  and  24-hour  forecasts  of  wave  heights  are  trans¬ 
mitted  to  Fleet  Weather  Central,  Rota,  Spain. _ 
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